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ABSTRACT 

Context. Transmission spectroscopy is a powerful technique for probing exoplanetary atmospheres. A successful ground-based ob¬ 
servational method uses a differential technique that uses high-dispersion spectroscopy, but it only preserves narrow features in 
transmission spectra. Broadband features, such as the remarkable Rayleigh-scattering slope from possible hazes in the atmosphere of 
HD 189733b as observed by the Hubble Space Telescope, cannot be probed in this way. 

Aims. Here we use the chromatic Rossiter-McLaughiin (RM) effect to measure the Rayleigh-scattering slope in the transmission 
spectrum of HD 189733b with the aim to show that it can be effectively used to measure broadband transmission features. The 
amplitude of the RM effects depends on the effective size of the planet, and in the case of an atmospheric contribution therefore 
depends on the observed wavelength. 

Methods. We analysed archival HARPS data of three transits of HD 189733b, covering a wavelength range of 400 to 700 nm. The 
radial velocity (RV) time-series were determined for white light and for six wavelength bins each 50 nm wide, using the cross¬ 
correlation profiles as provided by the HARPS data reduction pipeline. The RM effect was first fitted to the white-light RV time series 
using the publicly available code AROME. The residuals to this best fit were subsequently subtracted from the RV time series of each 
wavelength bin, after which they were also fitted using the same code, leaving only the effective planet radius to vary. 

Results. We measured the slope in the transmission spectrum of HD 189733b at a 2.5cr significance. Assuming it is due to Rayleigh 
scattering and not caused by stellar activity, it would correspond to an atmospheric temperature, as set by the scale height, of T = 
2300 + 900K, well in line with previously obtained results. 

Conclusions. Ground-based high-dispersion spectral observations can be used to probe broad-band features in the transmission spec¬ 
tra of extrasolar planets, such as the optical Rayleigh-scattering slope of HD 189733b, by using the chromatic Rossiter-McLaughiin 
effect. The precision achieved with HARPS per transit is about an order of magnitude lower than that with STIS on the Hubble Space 
Telescope. This method will be particularly interesting in conjunction with the new echelle spectrograph ESPRESSO, which currently 
is under construction for ESO’s Very Large Telescope, which will provide a gain in signal-to-noise ratio of about a factor 4 compared 
to HARPS. This will be of great value because of the limited and uncertain future of the Hubble Space Telescope and because the 
future James Webb Space Telescope will not cover this wavelength regime. 

Key words. Planets and satellites: atmospheres, individual: HD 189733b - Techniques: radial velocities - Methods: observational 


1. Introduction 

Transmission spectroscopy has proven to be a powerful tech¬ 
nique for probing the atmospheres of extrasolar planets ( |Char- 
bonneau et al.||2002t|Deming et al.||2013| ). Ground-based trans¬ 
mission spectroscopy is very challenging, even though tele¬ 
scopes with much larger collecting areas can be used than from 
space. Observations from the ground are hampered by the fact 
that a target is always seen through a varying amount of at¬ 
mosphere. Moreover, atmospheric circumstances are never per¬ 
fectly stable over a timescale of a few hours. In addition, in¬ 
strument stability is crucial, since ground-based observations 
are subject to changes in the gravity vector and/or field rota¬ 
tion. One way to mitigate these issues is by using ground-based 
high-dispersion spectroscopy. Eor example, detections of sodium 
and potassium have been presented in the optical(|Redfield et al. 
[20081 [Snellen et al.||2008l[Sing et al.||20nH|Wood et al.||2011t 


Jensen et al.||2QlH [Zhou & Bayliss|[2Q12t jWyttenbach et al. 
2015j ), and detections of carbon monoxide and water in the in¬ 


frared (Brogi et al. 


MT2l|Rodler et al.|2012t[de Kok et al.|2013t 


Birkby et al.|2013[ Lockwood et al.|2()14[|Brogi et al.|2014|). 


Unfortunately, the high-dispersion spectroscopic technique 
is insensitive to broadband absorption features since only rel¬ 
atively narrow features can be probed because all large-scale 
structures in the spectra are removed during the data analysis to 
account for the telluric variability. Therefore, broadband features 
can only be probed from the ground using multi-object spec¬ 
troscopy, which simultaneously observes the target and a num¬ 
ber of reference stars to perform differential spectrophoto metry 
to correct for atmospheric effects ( [Bean et al.|2010||2011| ). 


However, the most precise observations have been acquired 
from space. A particularly interesting broadband feature is the 
remarkable slope in the transmission spectrum of HD 189733b, 
as measured by [Pont et aL] ( |2008| ) and [Sing et aL] ( |2011b| ) us- 
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ing the Hubble Space Telescope (HST). [Pont et al.| ( r2013| ), who 
discussed the UV and optical spectrum based on a total of eight 
transits observed with STIS and ACS on the HST, found the ra¬ 
dius ofthe^lanetto_decreaseb^ .5% be twe en 345 nm and 775 
nm. [Lecavelier Des Etangs et al.lpOOS] ) and [Pont et al.| ( [MT3] ) 
interpreted this slope as caused by Rayleigh scattering by hazes 
at a temperature of ~ 1300 K, possibly increasing to ~ 2000 K 
at the shortest wavelengths probing the lowest pressure region 
in the atmosphere. It would be extremely valuable if there were 
another way to probe this feature from the ground. 

|Snellen| ( |2004| ) devised a method for probing the varying 
planet size as function of wavelength of a transiting exoplanet 
in a different way. The technique makes use of the Rossiter- 
McLaughlin (RM) effect. If the orbital plane of the planet is 
aligned with the spin of the star, the transiting exoplanet will first 
block light from the approaching part and then from the receding 
part of the stellar surface because the host star rotates. This ef¬ 
fect results in a wobble in the radial velocity (RV)of the star, first 
seen by [Rossit^ ( | 1 924| ) for eclipsing binaries, and recently ob- 
served for transiting planets in many systems (e.g. |Queloz et al. 


|2000[ ). To first order, the overall amplitude of this effect. Arm, is 
proportional to 


Arm ^ y sin i 


Rl 


DZ _ D. 

^star y 


( 1 ) 


where V sin i is the projected rotation velocit y of the star, 
is the planet radius, and Rstar is the stellar radius ( Haswell|2010| ). 
As can be seen, this amplitude depends on the effective radius of 
the planet, and since this (in case of observable transmission fea¬ 
tures) is wavelength dependent, subsequently the amplitude of 
the RM effect is also wavelength dependent and can be measured 
( Snellen|2004| ). The advantage of this method over conventional 
transmission spectroscopy measurements is that instead of us¬ 
ing the intensity of off-transit spectra as a reference for in-transit 
spectra, it depends on the changes in the profiles of the stellar 
spectral lines in the same on-transit spectra. For ground-based 
observation this makes the technique less likely to be influenced 
by Earth atmospheric effects, which turned out to be challenging 
for the conventional approach, because line shapes are less influ¬ 
enced by telluric absorption. Numerical simulations have been 
performed by preizler et al.|(|2009 ) to show the feasibility of this 
technique. Recently, [Czesla et al.| ( 2Q12| ) used a technique based 
on very similar principles to study the chromosphere of active 
planet host-stars. 

While [Snellen| ( [2004] ) advocated this technique to probe nar¬ 
row absorption features, such as those from sodium, we now 
realize that this method can be very powerful in probing wide 
broadband features, which are very challenging to probe from 
the ground. In this paper we test this technique on the exoplanet 
HD 189733b using archival HARPS data. In Sect. 2 the archival 
data and initial data analysis are described, and in Sect. 3 the 
models used to fit the RM effect are presented. In Sects. 4 and 5 
we present the fit to the data and discuss the final results. 


2. HARPS data and radial velocity extraction 

HD 189733b is a hot Jupiter with a mass of 1.15 Mj^piter that 
orbits a K dwarf in 2.2 days. The system parameters used for 
our analysis are presented in Table [T] It is one of the two most 
frequently studied and observed exoplanets, together with HD 
209458b. Its atmosphere has been observed both in trasmission 
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and during secondary eclipses (e.g.[Deming et al.|2006[[Huitson 
|et al.|201^|Wyttenbach et al.|2015|). 

We used HARPS archival data on HD 189733b, obtained 
under the programs 072.C-0488(E), 079.C-0828(A), and 079.C- 
0127(A). We chose not to use the data from 2006-07-29 because 
of the partial coverage of the transit (first half only) and because 
of the incomplete observation of the baseline. Instead, we used 
observations obtained during a total of three transits, w hich is the 
same dataset as was used by |Wyttenbach et ah ( 2015|), see t heir 
Table 1 for the observation logs. HARPS ( [Mayor et al.|2003] ) is a 
fibre-fed, cross-dispersed high-resolution (R ~ 115000) echelle 
spectrograph at the ESO 3.6m telescope, which through its sta¬ 
bility is ideal for measuring radial velocities. It covers wave¬ 
lengths from 378 nm to 691 nm over 72 orders. The data were 
taken during the night of September 7, 2006, consisting of 28 ex¬ 
posures of 600 seconds each, and on the nights of July 19, 2007 
and August 28, 2007, consisting of 38 and 40 exposures of 300 
seconds each, respectively. We retrieved the cross-correlation 
functions (CCFs) from the ESO archive. These CCFs are calcu¬ 
lated by correlating the data with a binary mask, both computed 
for each order of the spectrograph and for the full wavelength 
range, which are standard products of the ESO archive. For com¬ 
parison, we also used the values computed for the radial velocity 
using the full wavelength range and the errors on these values. 
We note that the data reduction as used by the ESO archive used 
different masks for different nights, that is, a K5 mask for the 
first night and a G2 mask for the others. We used the data for 
the first night reduced with the G2 mask with the HARPS Data 
Reduction Software by |Wyttenbach et ar|pQ15| ). 


2.1. Radial velocity extraction from the orders 

To determine the radial velocity as function of wavelength for 
each observed spectrum, we used the CCFs determined for each 
order of the HARPS spectrograph. We first combined the CCFs 
by averaging them, inside wavelength passbands of 50 nm cho- 
sen to b e the same as used by |Pont et al.| ( |2008] ) and |Sing et al.| 
( |2011b| ). These are shown in Table In this way, we can keep 
the same limb-darkening coefficients as used in their studies and 
compare them directly. We note that the fourth and fifth passband 
overlap by half in wavelength because this is where the individ¬ 
ual studies of [Pont et~ar] ( |2008| ) and |Sing et al.| ( |2QlTb| ) overlap. 
The second column of Tablej^shows the range of HARPS orders 
in each passband. There is no exact match with the passbands 
used by Pont et ah] ( [2008| ) and [Sing et ah] ( |2011b| ), but differ¬ 
ences are negligible. The orders were assigned to the different 
passbands based only on their central wavelength. HARPS order 
97 (central wavelength A = 631.06nm) was removed from our 
analysis because of strong telluric contamination from molecu¬ 
lar oxygen absorption. 

In this way, we obtained seven CCFs per spectrum for the 
seven different passbands. The radial velocity extraction was per¬ 
formed as in the HARPS Data Reduction Software, that is, by 
fitting Gaussians to the CCFs. The centre of the Gaussian is the 
RV of that passband. This approach of averaging the CCFs in a 
passband before determining the radial velocity is similar to the 
standard HARPS data reduction procedure, in which the CCFs 
from the different orders are first combined before they are fit 
with a Gaussian. 

The photon-noise-limited uncertainty of a radial ve locity 
measurement is discussed by [Hatzes & Cochran] ( |1992| ), who 
showed that 
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where ctrv is the uncertainty in the radial velocity, S is the flux, 
and R is the spectral resolution. The variable AT in practice is 
proportional to the number of spectral lines and their depths over 
the observed wavelength range. Since this is challenging to de¬ 
rive from first principles for the different passbands, we used the 
dispersion in the radial velocity time-series to scale the errors 
afterwards. 


3. Models for the Rossiter-McLaughlin effect 


We used the publicly available code AROME ( |Boue et al.|2013| ) 
to analytically compute the RM radial velocity anomaly as mea¬ 
sured with the CCF method. This code, written in C, was given 
the orbital planetary and stellar parameters (the projected rota¬ 
tion velocity Vsin(i), the semi-major axis a, the radius of the star 
Rstar, expressed in solar radii units, the inclination angle /, the 
mutual incl ination angle T, and the planet orbital period P, all 
taken from |Triaud et al.| ( |2009| ), see Table [TJ, and the times of 
observation, which it subsequently used to calculate the CCFs 
and the expected RV measurements from a Gaussian fit to these. 
The fraction of flux blocked by a planet crossing the surface of 
its host star depends on the limb darkening of the star. For the 
calculation with AROME we chose the non-linear law of IClaretl 
( |2000|) to matc h those used in the analysis of [Pont et~ar] ( 2008 1 
and jSing et al.| ( |2Qirb| ), that is. 


I(j^) 

/(I) 
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(3) 


For the white-light RV time seri es we used the qua dratic 
limb-darkening law coefficients as in |Triaud et al.| ( [20n9) ), who 
studied the RM effect for this planet using the same data set. All 
the limb-darkening coefficients used in our analysis are given in 
Table Since we are interested in the dependence of the RM 
effect on the planet effective size, for each passband we selected 
the limb-darkening law corresponding to those wavelengths and 
calculated the anomaly for 5000 values of the planet-star radius 
ratio (Rp/Rstar), between 0.105 and 0.205 with a step of 2 • 10“^. 
These boundaries were chosen arbitrarily with the purpose of in¬ 
cluding the measured planet-star radius ratio for this planet using 
the same HARPS data as |Triaud et al.| ( |2009| ) (0.1581 + 0.0005). 
All other parameters were kept fixed. We note that for the analy¬ 
sis we are only interested in the relative change of the Rp/Rstar as 
function of wavelength, not in absolute values. Therefore it is not 
relevant to fit for parameters such as the projected rotation veloc¬ 
ity Vsini, the mutual inclination angle T, the macroturbulence z, 
and the scaled semi-major axis <2/Rstar separately for each pass- 
band. These should be independent of wavelength and therefore 
can only result in a small vertical offset in the transmission spec¬ 
trum, not in a change of slope. 


4. Data analysis 

The aim of our data analysis is to determine the change in the 
effective radius of the planet over the seven wavelength pass- 
bands. The first step in the analysis was to remove the RV varia¬ 
tions due to the reflex motion of the host star around the planet- 
star barycenter. In addition, common mode errors that are ex¬ 
pected to be similar over the whole wavelength range covered by 
HARPS were removed, since they will not influence the change 
in Rp/Rstar as function of wavelength. The final step is the re¬ 
trieval of Rp/Rstar for the seven passbands. 


Table 1. Planet orbital and stellar parameters used to c alculate the 
Rossiter-McLaughlin (RM) effect with the AROME code ( |Boue et si] 
2013|l. The orbital planetary and stellar parameters are taken from and 
Triaud et ^ ( [2009 1 ). The remaining parameters, i.e. the line width of the 
non-rotating star j3o, the sigma of the Gaussian fit to the CCE (Tq, and 
the macro-turbulence parameter z, were estimated by us and only have 
a marginal effect on the resulting RM anomaly. 


Parameter 

value 

units 

reference 

a 

0.03120 

au 

Triaud et al. 

( 

2009) 

-^star 

0.766 

Re 

Ifiaud et al. 

( 

20091 

i 

85.508 

deg. 

Ifiaud et al. 

( 

2009 ( 

A 

-0.85 

deg. 

Ifiaud et al. 

( 

2009 ( 

Vsin(0 

3.10 

km/s 

Ifiaud et al. 

( 

2009 ( 

Po 

1.3 

km/s 

- 


cro 

3.3 

km/s 

- 


z 

4.0 

km/s 

- 


p 

2.218573 

days 

Triaud et al. 

(2009) 


Table 3. Eitted values of the systemic velocity Vsys and semi-amplitude 
of the reflex RV motion of the star, K for the three nights. 


Observation 

t^sys 

K 

Night 1 

-2.2431 +0.0002 

-0.202 ± 0.003 

Night 2 

-2.2366 ± 0.0008 

-0.190 + 0.007 

Night 3 

-2.2238 ± 0.004 

-0.205 + 0.005 


4.1. Reflex motion of the star 


Assuming that the orbit of HD 189733b is circular, the RV mo¬ 
tion of the host star Vstar can be modelled as a sinusoid with semi¬ 
amplitude K and period P given in Table plus a systemic ra¬ 
dial velocity Vsys- The exact values of Vsys differ from night to 
night because of the different stellar templates used to compute 
the CCFs (see Sect. 2). In addition, the RVs computed for the 
different passbands will also have small offsets because of the 
uncertainty in the reference RV for each order, see below. Stel¬ 
lar activity may also influence Vsys by spots covering parts of the 
rotating stellar disk. In addition, star spots rotating in and out 
of view may also result in small differences in K between the 
nights ( [Albrecht et al.||MT^ . We performed a least-squares fit 
of the out-of-transit RV data to obtain K and v^s for each night. 
These values are shown in Tableand in Fig.[2 


The reflex motion was subsequently removed from the 
white-light RV time series and from that of each passband. Small 
residual velocity offsets at a level of a few hundred meters per 
second were still present in the different passband RV time-series 
(since different parts of the templates were used). These offsets 
were also fitted and removed. 


4.2. Removal of common mode errors from the data 

Common mode errors that are identical across all passbands be¬ 
cause of residual instrumental effects, star spots (to first or der), 
stellar differential rotation, or other systematic effects (e.g. 1^ 
jaud et al.]|2009j ), can be removed from the data without harm. 
Therefore, we first performed a least-squares fit of the AROME 
models (Sect. 3) to the white-light RV time-series to determine 
the RM anomaly. This best fit was then removed from the white- 
light data to give the common mode errors on the time series, 
which were subsequently removed from the RV time series of 
each passband. To fit the white-light RV time series, a quadratic 
limb-darkening law was used with coefficients from jTriaud et al.j 


Article number, page 3 of[^ 


































































A&A proofs: manuscript no. hdl89_revised 


Table 2. Wavelength passbands, the corresponding HARPS orders, and limb-darke ning coefficients as used in our analysis. The first five pass bands 
and their non-linear limb-darkening coefficients were chosen to match those from Sing et ar](|2011b| ), while the last three match those from |Pont| 
|et al.| ( [2008] ). The coefficients for the entire HARPS range were taken from |Triaud et al.|p009| ) for a quadratic limb-darkening law. We note that 
order 97 (central wavelength A = 631.06 nm) was strongly affected by telluric mo' 
analysis. 


ecular oxygen absorption and was therefore removed from the 


Passband (nm) 

Orders 

cl 

c2 

c3 

c4 

370-420 

146-161 

0.5836 

-0.8102 

1.7148 

-0.5389 

420-470 

131 - 145 

0.5089 

-0.4084 

1.3634 

-0.5302 

470-520 

118-130 

0.5282 

-0.3141 

1.1931 

-0.4994 

520-570 

108-117 

0.6158 

-0.3460 

1.0695 

-0.4578 

550-600 

104-111 

0.4621 

-0.2003 

0.9450 

-0.4045 

600-650 

95 -96,98- 102 

0.5148 

-0.2774 

0.9429 

-0.4033 

650-700 

90-93 

0.5620 

-0.3531 

0.9524 

-0.4096 

Triaud 

90-161 

0.6355 

0.1488 





Orbital phase 


Fig. 1. White-light RV time series, plotted as blue dots for the first night, red squares for the second, and green triangles for the third night. The 
solid lines, with the same colours as used above, are the fitted reflex motion (K) and Vsys for each night. These values are shown in Table[^ 


(2009), who analysed the same HARPS data as here. These co¬ 
efficients are shown in Tabled 

To determine the best value of the radius ratio Rp/Rstar for the 
white-light RV time series, we computed the chi-squared;^^ for 
each model to the data, with the uncertainties rescaled such that 
the reduced of the best-fit model is equal to unity. The Icr- 
confidence interval was then determined by selecting the values 
of the radius ratio at which Ax^ = 1. 


The best fit, shown in Fig.[^ yielded a planet/star radius ratio 
of Rp/Rstar = 0. 1585 + 0.0007, in agreement with the value found 
by|Triaud et al.|(|2009|) for this planet based on the same HARPS 
data (0.1581 +0.0005). 
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Fig. 3. Fit of the Rossiter-McLaughlin effect for all the individual wavelength bins (from top left to bottom right). In the top panels of each figure, 
the RVs are shown in blue, computed for each bin and in red the fit to them; in the bottom panels of each figure the residuals of the fit are shown. 


4.3. Fit of the Rossiter-McLaughlin curve 

The resulting time series were fitted using the grid of models de¬ 
scribed in Sect. by minimization. The uncertainties on the 
measurements were rescaled such that for the best fitTreduced “ ^ ’ 
and the associated Icr errors were determined from the condition 
that = 1. We note that even if the limb-darkening coeffi¬ 
cients chosen for each passband were wrong by a quarter of a 
passband width, this would only contribute to the uncertainty on 
the retrieved /?p//?star value on a level of less than O.Scr. 


5. Results 

The least-squares fits to the RM effect in the individual pass- 
bands are shown in Fig. The scatter of the residuals ranges 
from 2 m/s for the 470-520 nm passband to 9 m/s for the 650- 
700 passband. The Rp/Rstar obtained for the different wavelength 
passbands are given in Table and are shown in Fig.[^ In Fig. 
I^we show the data of the fifth passband with the best-fit model, 
and that for 5cr smaller and larger planet/star radius ratios to il¬ 
lustrate the sensitivity of our fitting method. 

Star spots, unocculted by the planet, change the effective 
size of the star as function of wavelength. Therefore, they can 
have an equivalent effect on Rp/Rstar as a change in the ef- 


Table 4. Rp/Rstar values obtained for the different wavelength pass- 
bands. The quoted uncertainties are for the relative values, not for ab¬ 
solute measurements. The larger error at the longest wavelength arises 
because there are only four spectral orders contributing to that passband. 
The errors are for a Icr confidence level. 


Passband (nm) 
370-420 
420-470 
470-520 
520-570 
550-600 
600-650 
650-700 


-^p /^star 

0.1626 + 0.0015 
0.1610 + 0.0009 
0.1606 + 0.0007 
0.1582 + 0.0011 
0.1591 +0.0015 
0.1595 + 0.0012 
0.1618 + 0.0035 


fective size of the planet, which is what we wish to measure 
with transmission spectroscopy. [Sing et aT] ( |2011b| ) corrected 
for unocculted star spots by assuming a spot temperature of 
T'spot = 4250 + 250 K, causing a stellar flux reduction of 1% 
at 600 nm. We performed that same analysis. However, because 
our uncertainties in Rp/Rstar are an order of magnitude larger, 
this spot correction contributes only at the 0.1 cr level and was 
therefore ignored. 
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Fig. 3. Continued... 


5.1. Rayleigh-scattering slope 


As in [Siiig et al.| ( |2011b| ), [Pont et aL] ( [2008| ), and [Pont et aL] 
( 2Q13| ), we interpreted the optical transmission spectrum of 
HD 189733b as due to Rayleigh scattering. The Rayleigh scat¬ 
ter cross-section can be written as cr = (Tq (T/Tq)^, with a = -4, 
and therefore the slope of the planet radius as a function of wave¬ 
length is given by |Lecavelier Des Etangs et al.|p008| ) 




(4) 


where p is the mean mass of the atmospheric particles, taken as 
2.3 times the mass of the proton, k is the Boltzmann constant, 
and g the surface gravity. From the slope of the fit to the Rp/Rstar 
ratio as a function of the natural logarithm of the wavelength, an 
estimate of the atmospheric temperature (at pressures where the 
scattering takes place) can be derived. 

We performed a least-squares fit for the data as shown in 
Fig. weighting down the data points of passbands 520- 570 
nm and 550 - 600 nm by a factor 0.75 because they overlap by 
half. The best fit is shown in Fig. and coresponds to a gradient 
of -0.0064 + 0.0026 ln[A]“^, which is significant at a 2.5cr level. 
Assuming Rstar = 0.7667?© ( [Triaud et al.|2009| ), this gradient cor¬ 
responds to an atmospheric temperature ofT = 2300 + 900K 
(errors at Icr confidence level). 
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We compared our results with those of |Pont et al. ( 2013|), 
who performed the combined analysis of the data in [Sing et al 


( 2011b| ) and |Pont et al.| ( [2008 ). This is shown in Fig.[^ We note 
that we shifted the data from Pont et al.] ( |2013| ) up by 0.0002 in 
^p/^star to match our values. Our analysis is not meant to give 
precise absolute values of 7?p/7?star (which is infiuenced by the 
choices in stellar parameters such as projected rotational velocity 
and macro turbulence), only the relative change as function of 
wavelength. The uncertainties in our measurements are typically 
an order of magnitude larger than those of |Pont et~ST] ( |2013[ ). 


Pont et al. ( [2013 ) inferred two different atmospheric temper¬ 


atures for two separate wavelength regimes, below and above 
550nm. Above 550nm they found a gradient that corresponds 
to a temperature of about 1300K, while at shorter wavelengths 
the slope steepened corresponding to a temperature of about 
2000 K. As can be seen in Fig. the uncertainty in our data 
is too large to fit multiple slopes. The uncertainty in temperature 
found in our analysis agrees with both temperatures as found by 
[Pont et ar]p013| ). The temperature in the upper atmosphere of 
HD 1 89733b has also been estimated usign the Nal doublet line 
Huitson et al.j ( |2012| ) used HST data and found a temper- 


cores. 


ature of 2800 + 400 K, while jWyttenbach et al. ( |2015| ) found a 
temperature of 2500 + 400 K, using the same HARPS data set. 
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Fig. 2. Top panel: the white-light RV time-series with the best fit to 
the Rossiter-McLaughlin anomaly. The fit yields a planet/star radius ysl- 
ti o of Rp/Rstar = 0.158 5 + 0.0007, in agreement with the measurement 
of |Triaud et al.| ( [200^ . In the bottom panel the residuals to the fit are 
shown. These residuals contain information about possible star spots 
crossed by the planet during a transit and possible effects due to stellar 
differential rotation or macroturbulence. These residuals were subse¬ 
quently removed from the radial velocity time-series of the individual 
passbands. 



Fig. 4. Upper panel: the radial velocity time-series for the fifth passband 
(550-600 nm) and the corresponding best-fit model (blue dots and red 
line, respectively). The lower panel shows the residuals to the fit, to¬ 
gether with the 5cr upper and lower limit of the Rp/Rstar ratio (solid red 
and green lines, respectively). 


5.2. Other interpretations of the slope 


In our analysis we assumed that the slope in the transmission 
spectrum of HD 189733 b is due to Rayleigh scattering, as was 
also ass umed by [Pont et al.|p008]),|Lecavelier Des Etangs et aL| 


( [2008| ), Sing et al.| ( |2011b| ), and |Pont et al.| ( |2013| ). However, in 
recent literature it has been advocated that star spots or plages 
might also result in the observed slope in the transmission spec¬ 
trum. [McCuIIoughet^r|(|20^ reinterpreted the available mea¬ 
surements with a clear planetary atmosphere. They found that 
an unocculted spot fraction of 4% can mimic the increase in ef¬ 
fective planet size towards shorter wavelengths. Occultation of 


Fig. 5. Resulting planet-star radius ratio Rp/Rstar as function of wave¬ 
length, indicating the Icr error bars. The solid line indicates the least- 
squares fit of a Rayleigh-scattering slope to the data, using Eq.|^ The 
dashed lines indicate the Icr error margin on the slope. Note that the 
abscissae axis is in a logarithmic scale. The linear fit yields a gradi¬ 
ent of -0.0064 + 0.0026 ln[A]"^ corresponding to a temperature of 
T = 2300 + 900K. 



Wavelength (A) 


Fig. 6. Comparison of our Rp/Rstar measurements (blue squares) with 
those of [Pont et~n^ ( |2013| ) (red circles). The Icr error bars are indicated. 
The solid and dashed lines are the same as in Fig.|^ 


stellar plages has also been claimed to be able to give rise to the 
slope in the spectrum of HD 189733b by |Qshagh et al.J ( [2014| ). 


We note that the chromatic RM effect as used in our analysis to 
determine the effective planet radius as function of wavelength 
is affected by star spots and/or plages in a similar fashion as clas¬ 
sical transmission spectroscopy. Therefore, this analysis cannot 
contribute to this discussion at this stage. 

However, it is worth noting that an increase in planetary ra¬ 
dius towards UV wavelengths has been detected for a number 
of planets, for instance, GJ 3470b ([Nascimbeni et al.|201 3} [Bid- 


Idle et al.|2Q14l), W ASP-6b ( [Nikolov et al.|2015| ), and ^\^P-31b 


( [Sing et al.|2015| ). Interestingly, the host stars of all these plan¬ 
ets and that of HD 189733b are active stars. This might indicate 
that unocculted star spots or plages may be the cause of this ef¬ 
fect. [Qshagh et al.| p014| ) also advocated this for the particular 
case of GJ3470b. A more detailed investigation is required to 
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6. Conclusion 
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Fig. 7. Simulated Rp/Rstar for five transits observed with the future 
ESPRESSO spectrograph on the VLT (blue squares). Observatio ns of 
five transits by the STIS spectrograph on HST as presented by |Pont| 
|et al.| ( [20T^ are indicated by red circles. These are shifted slightly red- 
ward for clarity. All error bars indicate Icr uncertainty intervals. These 
are larger by between 1.5 and 2 times for the ESPRESSO than for the 
STIS observations. 


We used the chromatic Rossiter-McLaughlin technique to probe 
the slope in optical transmission spectrum of HD 189733b us¬ 
ing archival HARPS data. We showed that this technique can 
be used to measure this slope and found a change in radius ra¬ 
tio Rp/Rstar of -0.0064 + 0.0026 ln[A]“^ (2.5cr). Assuming this 
is due to Rayleigh scattering, it corresponds to an atmospheric 
temperature ofT = 2300 + 900K. 

The precision achieved with HARPS per transit is about an 
order of magnitude lower than that with STIS on the Hubble 
Space Telescope. This method will be particularly interesting 
in conjunction with the new echelle spectrograph ESPRESSO, 
which currently is under construction for ESO’s Very Large 
Telescope and will provide a gain in signal-to-noise ratio of 
about a factor 4 compared to HARPS. This will be of great value 
because of the limited and uncertain future of the Hubble Space 
Telescope and because the future James Webb Space Telescope 
will not cover this wavelength regime. 
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understand the relationship between stellar activity and slopes in 
transmission spectra of orbiting planets. 


5.3. Future observations with ESPRESSO 


We showed that with three transit observations of HD 189733b 
with HARPS, we obtained a relative precision of the planet ra¬ 
dius as a function of wavelength, which is about an order of 
magnitude lower than that a chieved for five transits with STIS 
at the HST ( [Pont et al.|2Q13| ). Since the chromatic RM technique 
is photon limited, a significant increase in precision can be ex¬ 
pected from the forthc oming ESPRESSO spectrograph at ESO’s 
Very Large Telescope ( |Pepe et al.|2Q13| ). 

We simulated future ESPRESSO observations by scaling the 
current accuracy achieved with HARPS by the increase in col¬ 
lecting area of the telescope and the increase in throughput of 
the ESPRESSO spectrograph, also taking into account the slight 
increase in spectral resolution to R = 120000. The photon- 
noise contribution to the radial velocity error is, according to 
Eq-i ctrv oc S where S is the received flux from the 

star, and R the resolving power. Eigure shows the results of 
our simulation for flve transits observed with ESPRESSO (blue 
squares) with their associated Icr unc ertainties. In red ci rcles are 
shown the data points as presented by [Pont et al.| ( |2013] ) for also 
flve transits observed with STIS. The expected uncertainties for 
ESPRESSO are about 1.5-2 times larger than for STIS. This is 
particularly important, taking into account the limited and uncer¬ 
tain future of the Hubble Space Telescope, and the fact that the 
future James Webb Space Telescope will not cover this wave¬ 
length regime. 

We do note that the amplitude of the RM effect depends on 
the spin of the star and its angle with the planet’s orbit. There¬ 
fore, not all transiting planetary systems, for instance those or¬ 
biting slow rotators, are equally favourable for this technique. 
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